High Magnetic Field NMR Studies of LiVGe 2 O e , a quasi 1-D Spin S = 1 System. 
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We report 7 Li pulsed NMR measurements in polycrystalline and single crystal samples of the quasi 
one-dimensional S — 1 antiferromagnet LiVGe2 06, whose AF transition temperature is Tn — 24.5 K. 
The field (Bo) and temperature (T) ranges covered were 9-44.5 T and 1.7-300 K respectively. The 
measurements included NMR spectra, the spin-lattice relaxation rate (TV -1 ), and the spin-phase 
relaxation rate (T 2 _1 ), often as a function of the orientation of the field relative to the crystal 
axes. The spectra indicate an AF magnetic structure consistent with that obtained from neutron 
diffraction measurements, but with the moments aligned parallel to the c-axis. The spectra also 
provide the T-dependence of the AF order parameter and show that the transition is either second 
order or weakly first order. Both the spectra and the T7" 1 data show that Bo has at most a small 
effect on the alignment of the AF moment. There is no spin-flop transition up to 44.5 T. These 
features indicate a very large magnetic anisotropy energy in LiVGe2C>6 with orbital degrees of 
freedom playing an important role. Below 8 K, T^ 1 varies substantially with the orientation of Bo 
in the plane perpendicular to the c-axis, suggesting a small energy gap for magnetic fluctuations 
that is very anisotropic. 



PACS numbers: 75.30.Kz, 75.50.Ee, 76.60.-k 

I. INTRODUCTION 

Recently, a new quasi 1-D spin 5 = 1 sys- 
tem, LiVGeaOg,, has been the object of intensive 
experimental™ 2 ™ and theoretical investigations.oBtl It 
has an antiferromagnetic phase transition at about 25 K 
and the expected Haldane gap is either absent or strongly 
suppressed. Quantum chemistry calculations!!] indicate 
that a second-order splitting Acf of the t2 g orbitals may 
play a dominant role in this system. Our new measure- 
ments indicate that. Acf might be much smaller than 
previously thought, El leading to a large uniaxial magnetic 
anisotropy and orbital fluctuations. 

It has been established by neutron diffraction mea- 
surements that the low temperature phase .has a rather 
simple, long-range antiferromagnetic order.H In this pa- 
per, we report the results of a number of different NMR 
measurements on this material. We also address several 
important questions about the phase transition which re- 
mained open previously, including the order of the phase 
transition, the size and origin of the energy gap in the 
magnetic excitation spectrum below the Neel tempera- 
ture, and the orientation of the magnetic moments in 
the antiferromagnetic phase. Many of the results re- 
ported here were obtained on powder samples. Some 
of the more recent measurements were made on single 
crystal samples. 

The 7 Li NMR measurements we report include NMR 



spectra, the spin-lattice relaxation rate T± , and the 
spin-spin relaxation rate T 2 _1 , at magnetic fields B be- 
tween 9.0 and 44.5 T and temperatures T over the range 
1.7 - 300 K. In spite of various attempts to observe the 
resonance signal of 51 V nuclei above the transition and 
at the lowest temperatures in the AF-phase, only a tiny 
spurious signal could be detected in the polycrystalline 
sample and no signal at all was found in the single crystal 
samples. The 9.0 T measurements were made at UCLA, 
the measurements between 23 and 44.5 T were done at 
the NHMFL in Tallahassee, and measurements at 12 T 
were performed at the GHMFL in Grenoble. We have 
extended previous NMR measurements on a polycrys- 
talline powder samplco to much lower temperatures as 
well as to much higher magnetic fields. Furthermore, we 
present the first NMR measurements on LiVGe206 single 
crystals as a function of the polar and azimuthal angles, 
which give new insides on the low temperature behavior 
of this system, where orbital degrees of freedom seem to 
play an important role. 

This paper is organized as follows. First, we describe 
the preparation of the samples and the measurement pro- 
cedures. Then, we present the experimental results and 
a partial interpretation of some of them. In the subse- 
quent discussion we address issues concerning the mag- 
netic structure, the phase transition, the relaxation rate 
and the influence of orbital degrees of freedom. 
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II. SAMPLES AND EXPERIMENTAL 
METHODS 

The LiVGe2C>6 powder sample was prepared as de- 
scribed by Millet et alE The single crystal samples were 
synthesized at the Centre d'Etudes Nucleaires in Greno- 
ble using a flux of Ge02:Li2B407 with the molar ratio 
8:1. After reducing the V2O5 with H2 and using a slow 
increase of the temperature up to 720 C, the compound 
LiVGc206 was obtained with a thermal treatment of the 
components up to 800 C under Ar with 2 % vol H 2 . Then, 
a mixture of 70 %wt of the flux and 30 %wt of the com- 
pound was put in platinum crucible and heated for 1 day 
at 970 C. After that, it was slowly cooled at the rate 
of 2 C per hour to 780 C, after which the power to the 
furnace was switched off and allowed to cool to room tem- 
perature. Finally the products were washed with boiling 
water. Pale green needles were obtained, the maximum 
size of which was approximately 1 mm x 0.10 mm x 
0.050 mm. The typical dimensions of the samples used 
for our NMR measurements were 700 /iin x 100 fim x 
50 /im, which corresponds to a mass of about 15 fig and 
~ 3 • 10 16 7 Li spins. The small NMR coils used for most 
of the single crystal work were a few turns of 25 /Ltm 
diameter insulated copper wire wound tightly onto the 
sample. 

LiVGe2 0g crystallizes in the monoclinic system, space 
group P2i/cn The chains of VOg octahedra are paral- 
lel to the c-direction and are connected to their neighbor 
chains only by two Ge04 tetrahedra. There is a very 
small coupling perpendicular to the chains. The Vana- 
dium atoms are located in distorted oxygen octahedra 
and the three i2g orbitals are split into a low-lying dou- 
blet (d xy , d yz ) and a single orbital (d xz ) at an energy Acf 
above the doublet. 

All of the NMR results reported here were performed 
on the 7 Li nuclei using standard spin-echo techniques car- 
ried out with a spectrometer and probes built at UCLA. 
The NMR spectra were obtained by frequency-shifted 
and summed Fourier transform processing*] with fixed 
applied magnetic fields between 9.0 and 44.5 Tesla. Ro- 
tation of the field alignment about one axis during the 
measurements was done by placing the sample and NMR 
coil on a goniometer platform whose orientation was con- 
trolled from the top of the probe. Further rotation about 
a second axis perpendicular to the goniometer rotation 
axis was carried out by changing the placement of the 
coil and sample on the goniometer platform when the 
probe was out of the cryostat. We estimate that the ab- 
solute accuracy of the corresponding angle settings was 
approximately ±10 deg and that the precision in chang- 
ing the angle with the goniometer was ±0.5 deg. Part of 
the uncertainty in the absolute angle was associated with 
the small size of the samples and part of it from thermal 
contraction in the goniometer control upon cooling from 
room temperature to low temperatures. 

An unsuccessful attempt was made to observe the 51 V 
NMR signal in a single crystal sample. A thorough search 



was done by sweeping the resonance frequency in the 
absence of any external field at 4.2 K as well as sweeping 
B (aligned along the c-axis) between and 14 T with T 
in the range 1.5 K-5 K at the fixed frequencies 200, 300 
and 550 MHz. We attribute the lack of a signal, which 
otherwise should have been rather intense, to values of 
T\ or T2 that were less than the approximately 2 \is dead 
time of the NMR spectrometer. It may be that extending 
such measurements to much lower T will reveal this V 
signal. 

The 7 Li T^ 1 measurements were performed by first 
rotating the nuclear magnetization out of equilibrium 
by a short saturation chain of rf-pulscs, then waiting a 
variable recovery time, t, and finally measuring the in- 
tegrated spin-echo intensity, m(t). As discussed below, 
the quadrupolar splitting of the NMR line is very small, 
about 15 kHz, so that a single exponential form is ex- 
pected for m(t) as long as all parts of the sample have 
the same value of T£~ . To monitor any deviation from 
the single exponential behavior, we used a stretched ex- 
ponential to fit our data: 

m(t) = moo + (mo - moo) cxp(-t/Ti) /3 , (1) 

where tuq and moo are the nuclear magnetization at 
t = just after the saturation sequence and the equi- 
librium magnetization, respectively. The fit parame- 
ter T\ is the single time that characterizes the recov- 
ery of the magnetization. It is the time for the quantity 
[toqo — m(t)] /[moo — mo] to decay to 1/e. The exponent (3 
reflects the width of the distribution of relaxation rates. 
For j3 = 1, it corresponds to a single exponential and as 
[3 decreases from 1, it represents a progressively broader 
distribution. The T^ 1 measurements were done at Lar- 
mor frequencies between 149 and 762 MHz and applied 
magnetic fields between 9.0 and 41.5 T. For some of the 
measurements at high magnetic fields below 3 K only the 
beginning of the recovery curves were measured and the 
parameter Toqo was set using values from the measure- 
ments between 10 and 3 K and the inverse temperature 
dependence of moo- 

Our r 2 _1 measurements were done at 148.981 MHz in a 
field of 9.0 T. The pulse sequence used was a ir/2 prepa- 
ration pulse applied to TOqo followed a time r later by a 
second pulse whose angle was set to maximize the am- 
plitude of the echo. The integral of the spin echo signal 
was recorded as a function of r. The decay of the signal 
was analyzed using the function: 

m(2r) = TO(0)exp-(2T/T 2 ) /3 , (2) 

where j3 is the stretched exponential parameter between 
1 (exponential decay) and 2 (gaussian decay). For the 
powder sample the [3 parameter was usually left free dur- 
ing the fit, and resulted in values around 1.4. The spin- 
echo amplitude for the single crystal measurements was 
modulated by the quadrupolar interaction, which caused 
strong deviations of m(2r) from an exponential decay. In 
this case, the value of (3 used in the analysis was fixed at 
2. 
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III. EXPERIMENTAL RESULTS 

A. NMR spectra of polycrystalline LiVGe20s 

Figure |] shows two NMR spectra of the powder sample 
in the paramagnetic regime, i.e. at T > Tn w 25 K. The 
experimental points are indicated by the symbols and the 
solid and dotted lines are fits to a model of a polycrys- 
talline powder in, the presence of an axially symmetric, 
anisotropic shift as discussed below. 
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FIG. 1: 7 Li NMR spectra in the paramagnetic phase of 
polycrystalline LiVGe2 0e. The dotted line is the expected 
position of 7 Li in a reference compound like LiCl. The solid 
and dashed lines are simulations (see text). 

For the simulation of the asymmetric NMR spectra 
in the paramagnetic regime we assumed each V ion to 
have a moment along the applied magnetic field, whose 
magnitude is independent of the orientation. Hence, the 
anisotropy of the g- factor was not taken into account. 
It was, however, verified that when 1 > gj_/g\\ > 0.5, 
results similar to those for an isotropic <?-factor are ob- 
tained. The corresponding magnetic field at the Li sites 
was then calculated by adding all the dipole contributions 
of the V ions in a sphere of about 5 nm diameter around 
the Li ion. It was verified that modifying the diameter 
of the sphere does not have any effects on the results. 
The hyperfine field at the Li sites cannot, however, be 
fully accounted for by assuming a purely dipolar field of 
the V moments. An additional isotropic hyperfine cou- 
pling of the order of 0.048 T//ib is needed. The latter 
may arise in a manner similar to the superexchange in- 
teraction. By assuming a randomly distributed powder, 
i.e., the direction of the applied magnetic field is point- 
ing along all possible directions of the unit sphere, the Li 
spectra were then simulated at different temperatures. 
From the simulations of the measured spectra, one ob- 
tains for the average component < M% > of the magnetic 
moments along the direction of the applied field 0.048 fiB 



at 300 K and 0.080 /ib at 53 K, respectively. Their ratio 
is what one expects torn the temperature dependence of 
the dc-susceptibility.H The simulations, shown in Fig. |l| 
are based on only three parameters: i) an isotropic hy- 
perfine coupling, which is the same for all the data, ii) 
the size of the magnetic moment on the V ions, whose 
temperature dependence follows the dc-susceptibility and 
iii) a gaussian broadening function. In view of the small 
number of parameters the fits agree fairly well with the 
measured data. Some deviations are observed in the low 
frequency part of the signal; their origin is not yet un- 
derstood. 
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FIG. 2: 7 Li NMR spectra of polycrystalline LiVGe 2 6 at 
9.0 T for T between 15 and 27 K. 

The 9.0 T NMR powder spectra near the transtion 
and in the ordered phase are shown in Fig. |^. They show 
a continuous transfer of spectral weight from a narrow 
line in the paramagnetic phase to a broad signal in the 
antiferromagnetic (AF) phase that occurs over a narrow 
temperature range. As seen in Fig. [|, which shows the 
fraction of the intensity in the AF phase, both phases 
coexist in a temperature range of about 1.5 K around 
24.4 K. 

The order parameter of the AF phase is the magnitude 
and polarization of the AF moments associated with the 
V atoms. They generate a corresponding magnetic field 
Bi at the i-th 7 Li sites, which can be calculated for a 
given AF moment configuration. Depending on the con- 
figuration of the AF state, Bj may have a sequence of 
values at different Li sites or be given by a single value 
for all of them. At the high magnetic fields used in our 
experiments, the 7 Li spins probe the static order parame- 
ter through the shift in their spectrum, which is given by 
the component of B; that is parallel to the applied field 
Bo- For a randomly-oriented polycrystalline sample, the 
NMR spectrum in the AF phase depends on the response 
of the AF polarization to the varying orientation of the 
external field. 

The NMR spectra in the AF phase below 23 K have a 
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FIG. 3: Fraction of the total NMR intensity in the polycrys- 
talline LiVGe206 spectra at 9.0 Tesla which is attributed to 
the antiferromagnetic phase. 

broad, nearly rectangular shape. This shape is expected 
for a randomly oriented powder spectrum if Bi has the 
same magnitude at all 7 Li nuclei, is parallel or antipar- 
allel to a single crystalline direction and maintains the 
same orientation with respect to the crystalline axes for 
all orientations of Bo- For Bi <C Bq, the field at the 
nuclei is given by Bq + Bi cos(0), where 9 is the polar an- 
gle in spherical coordinates. Then, the frequency shift is: 
/ = r yBiCos(6). The probability (dN) that a particular 
value of 9 occurs is: dN — l/2sm(9)d9 and the density 
of states in the powder pattern is: 
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27-Bi 



(3) 



i.e., constant and results in a rectangular shape for the 
spectrum. 

Since this shape is indeed observed, the powder spectra 
indicate that the direction of the internal magnetization 
is not affected by the orientation of the external magnetic 
field. From these measurements (see Fig. |J) we obtain the 
value Bi = 0.106 T. 

In Fig. | the part of the NMR linewidth (Ai/(T)) pro- 
portional to the AF order parameter is plotted as a func- 
tion of temperature. The contribution of the width in the 
paramagnetic phase is subtracted in quadrature from the 
AF contribution using 

Av = (Av\ F - A^/2 , ( 4 ) 

where Avaf and Az/p are the HWHM of the spectra in 
the AF phase and in the paramagetic phase at 26 K, 
respectively. The solid line corresponds to a power law 
behavior Au(T) oc (T N - T)°- 4±0 05 . The onset of the 
broadening of the spectra occurs at 25.04 K, which we 
identify as the Neel temperature Tn for the polycrys- 
talline sample. 
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FIG. 4: Frequency shift Av from the NMR spectra of poly- 
crystalline LiVGe2 06- The width in the paramagnetic phase 
at 26 K has been subtracted. The solid line is a fit to the 
data (see text). 



B. NMR spectra of LiVGe2 0s single crystals 

The coordinates shown in Fig. |s| will be used to discuss 
our measurements on the single crystal samples. They in- 
clude the crystalline axes a, 6, and c, the cartesian axes 
x, y, and z, and the spherical coordinates 9 (polar an- 
gle) and (azimuthal angle). X-ray measurements^ have 
shown that the long dimension of LiVGe20g single crys- 
tals is along the crystallographic c-direction. For all of 
our single crystal NMR measurements, the value of <fi is 
close to zero, unless specified otherwise. 
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FIG. 5: Definition of various coordinates for a typical 
LiVGe2 06 single crystal. 

Figure ^ shows NMR spectra of a LiVGe2C>6 single 
crystal at temperatures between 23 and 26 K for 9 w 60°. 
The NMR spectrum is a single line in the paramagnetic 
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phase and two lines in the AF phase that correspond to 
the two magnetically inequivalent 7 Li sites. As for the 
powder spectra, both phases coexist over a small range 
of T. However, this range is substantially narrower with 
a value ~ 0.3 K for the single crystal. 
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FIG. 6: Single crystal 7 Li NMR spectra of LiVGe 2 6 near 
Tn. The needle direction, which corresponds to the crystallo- 
graphic c-direction, is aligned 60° from the applied magnetic 
field. 

As for the powder sample, the single crystal linewidths 
are rather broad near Tn. Just above the transition, 
Ar/(T) is about 36 kHz compared to 19 kHz at 38 K. Since 
the dc-susceptibility increases with increasing T in this 
range, the opposite change in the linewidth may indicate 
that there are slow, short range fluctuations above the 
phase transition that enhance the linewidth somewhat. 
Furthermore, in the whole considered temperature range, 
the linewidths of the single crystal signals are consistent 
with the width of the broadening function used to fit the 
spectra of the polycrystalline sample. The large NMR 
line width seen just below Tn may also be caused by a 
distribution of Tn in the sample. 

In Fig. 0, Az/(T), representing the splitting of the two 
peaks in the 7 Li NMR spectrum, is plotted as function 
of T for B a = 9.0, 41.5 and 44.5 T, with 9 = 30° at 9.0 T 
and approximately 20° for the measurements at 41.5 and 
44.5 T. The data obtained at 41.5 and 44.5 T have also 
been multiplied by 1.15 and 1.28, respectively, to bring 
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FIG. 7: Frequency shift Au of the Li peaks in the AF phase 
of LiVGe2 0e from the peak in the paramagnetic phase as a 
function of T. For comparison, the measurements at 41.5 and 
44.5 Tesla have been multiplied by the factors 1.15 and 1.28 
respectively. The solid line is a fit to the data (see text). 



the curves together. This indicates a somewhat reduced 
value of the order parameter for the measurements in 
higher fields, an aspect that will be discussed later. For 
the same reason, the values of T have been reduced by 
0.6 K and 1.6 K for the measurements at 41.5 T and 
44.5 T, respectively. We attribute these adjustments in 
T as being caused by differences in the instrumentation, 
which resulted in an uncertainty of about 1 K in Tn for 
the measurements at these higher values of Bq. 

Because of the large values of Bq, a substantial reduc- 
tion of Tn proportional to Bq is expectedil^ However, 
as can be seen in Fig. |7[ no such reduction is observed. 
This seems rather surprising because <?/xb-B/£;b = 54 K, 
for B = 44.5 T with g — 1.79, i.e., about twice Tn- It may 
be, however, that the value of g — 1.79, obtained from the 
paramagnetic susceptibilityEl is too large because, as wiU 
be discussed later, the measurements of Lumsden et al.,Q 
as well as our own measurements, show that the value of 
the ordered moment is only about 1.14 ^b- Furthermore, 
it is evident that a large change in Bq has no significant 
effect on the T dependence of the order parameter in the 
ordered phase, as the three curves coincide. The solid line 
is a power law fit to the data just below the transition 
with Az/(T) cx (T N - T)°' 3±0 05 . The value of the expo- 
nent is somewhat smaller than that 0.4 obtained for the 
polycrystalline sample. Qualitatively, such a reduction 
is consistent with the smaller distribution of Tn inferred 
from the smaller temperature range over which a coex- 
istance of two phases is indicated in the single crystal 
samples. 
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C. Tf 1 (T, B) for polycrystalline LiVGe 2 O e 



In this section we describe our spin-lattice relaxation 
rate measurements at 9.0 and 23 Tesla performed on the 
polycrystalline sample. Below Tn the measurements were 
done at the same frequency and magnetic field as above 
the transition. But because of the very broad linewidth 
of about 3.5 MHz in the AF-phase, only a small part of 
the NMR spectrum was covered by the RF-pulses. The 
nuclei near resonance were those with ks 90°; i.e., Bo 
close to lying in the azimuthal plane. The angle <f> is 
distributed randomly over 2tt. 
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FIG. 8: Recovery of the magnetization during a T x mea- 
surement in polycrystalline LiVGe2 0e. The solid line is the 
fit of Eq. [j] for (3 = 0.92. Inset: Plot of (3 as function of T at 
9.0 and 23 T. 



Figure g shows a typical magnetization recovery curve 
(filled circles). The best fit using Eq. [I] (solid line), which 
yields (3 — 0.92, is an excellent fit to the data. The inset 
of Fig. H shows that the T-dependence of j3 is essentially 
the same at both 9.0 T and 23.0 T. Except for T close 
to Tn, above 12 K the values obtained for (3 are close to 
1 and correspond to a relatively narrow distribution of 
Tf . Near the transition and below 12 K, (3 considerably 
deviates from 1, which indicates a substantially broader 
distribution of Tf 1 at these temperatures. We attribute 
the broadening of the distribution of Tf 1 close to Tn to 
the distribution of Tn for the different parts of the sam- 
ple mentioned earlier. As discussed below for the single 
crystal measurements, the large deviation from (3=1 
seen below 12 K reflects a large, unexpected dependence 
of Tf 1 upon (j). 

The Tf 1 values obtained for 9.0 T and 23.0 T are 
shown in Fig. ^. Where (3 is substantially less than 1, 
a wide distribution of values for Tf 1 is present. There- 
fore, the plotted value is the one that represents the single 
recovery rate that characterizes this distribution. Never- 
theless, from its T-dependence useful information on the 
dynamics of our system can be obtained, even at the low- 



10- r 



t-r io° j 



10- 1 



10" z ! 



10° I 



LiVGe -,0 6 polycrystal 



O 9.0 Tesla 
• 23 Tesla 



o 



A* 



9 



150 

is 100 

50 



20 22 24 26 28 

rot) 



uv io 2 

temperature T (K) 



FIG. 9: 7 Li Tf 1 as function of T in polycrystalline LiVGe 2 6 
at 9.0 (open circles) and 23 Tesla (filled circles). 



est values of T. For the 9.0 T measurements, well above 
Tn, Tf 1 depends only weakly on T, slowly increasing by 
about a factor 2 between 40 and 200 K and remaining 
almost constant between 200 and 300 K. Therefore, as 
previously reportedBcl, there is no indication of a Hal- 
dane gap in this quasi 1-D system. Below 40 K, Tf 1 
increases rapidly to a maximum value of about 140 s _1 
at Tn, presumably due to critical fluctuations near the 
transition. 

Below 23 K, Tf 1 drops very rapidly by about 6 decades 
to a value near 3 • 10~ 4 s _1 at 1.8 K. In the low T regime, 
the effect of B on Tf 1 (T) is very weak. Because the 
value of g^sB/k-B for 23 T is close to Tn this result was 
unexpected, as was the small influence of B$ on the spec- 
tra shown in Sec. Ill B] which are also affected only weakly 
by B . 



D. Tf X (T) of a LiVGe 2 6 single crystal 



Figure [l(j shows measurements of T l as a function of 
T for a single crystal of LiVGe 2 6 for = 90°, 60°, and 
30°, with (j> ~ 0°. Over the whole temperature range, 
no significant deviations from f3 = 1 were observed. This 
behavior indicates that unlike the polycrystalline sample, 
there is no distribution in Tf i.e., the relaxation follows 
a single exponential. In principle, this result should make 
a detailed interpretation of the results more straightfor- 
ward. 

The general behavior of Tf 1 (T) is similar to that of 
the polycrystalline sample. From 200 K to 40 K the 
relaxation rate slightly decreases by about a factor of 
two, but subsequently increases by more than one order 
of magnitude towards Tn, and decreases rapidly below 
T N . The inset of Fig. [lC] shows a log-log plot of Tf 1 
for « 90° close to T N as a function of |T - T N | with 
Tn = 24.45 K. Close to the transition, the data both 
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FIG. 10: 7 Li rf 1 as function of T in a single crystal of 
LiVGe2C>6. The three values of the polar angle 9 are 90° 
(open circles), 60° (open triangles), and 30° (filled circles). 
Inset: Tf 1 as function of \T — 24.45| above (open circles) and 
below (filled circles) T N for 9 « 90°. 



below and above Tn, fall on the same curve, given by 

Tf 1 oc |T- 24.45|~°- 55 . (5) 

Although the exponent —0.55 is expected to reflect the 
critical behavior of the AF transition in this material, it 
should be interpreted with caution. Because the width 
of the transition shown by the coexistence of both phases 
(Fig. ^|) over a range of 0.3 K may indicate a distribution 
of Tn in the sample, it may be that the exponent is really 
a lower limit on the rate of the divergence on approaching 



on T has a crossover from 9 being approximately inde- 
pendent of 9 below 8 K to a strong dependence above 
10 K. This behavior is shown in more detail in Fig. O, 
where Tf 1 is plotted as function of 9 for several values of 
T both above and below 10 K. Above 10 K, Tf 1 is well 
described by 



Tf 1 ^) =A(T) [35 sin 2 (9) + 5] 



(6) 



where A determines the mag nitude of Tf \ The mea- 
surements at 5 and 3 K do not show any dependence on 
9. At 2.2 K, the moderate ^-dependence of Tf is not 
well enough established to draw useful conclusions. 
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FIG. 12: 7 Li Tf 1 as function of 1/T in a LiVGe 2 6 single 
crystal at 28.5 T and the polar angle 9 = 90° for two az- 
imuthal angles <f>\ = —45° (filled circles) and <f>2 = 0° (filled 
squares) and at 9.0 T at 9 = 90° (open circles). The solid 
and dashed lines indicate energy gaps of 2.6 and 14.7 K re- 
spectively. The dotted line is a guide to the eye. 
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FIG. 11: Angular dependence of Tf 1 in a LiVGe206 single 
crystal at five temperatures between 35 and 2.2 K. The solid 
line is a fit to the data at 35 and 13 K (see text). 

Figure |l^ also shows that the 6*-dependence of T 1 _1 (T) 



Now we turn to measurements in which (f> was varied 
and 9 is held fixed at 90°; i.e., Bo was rotated in the 
azimuthal plane. For these measurements, the direction 
corresponding to 4> — is always the same, but the lo- 
cation of 4> = in the azimuthal plane is not known. 
Figure |l| shows Tf 1 at 6 = 90° as a function of 1/T for 
cf>i = —45° and rfe = 0° at 28.5 T and, for comparison 
at 9.0 T and cf> = 0°. In addition the cf> dependence of 
Tf 1 for 4.5 and 1.7 K at 6 = 90° and 28.5 T is shown in 
Fig. [l3|. A huge change of two decades in Tf is seen for 
a variation of only 30° in <j>. 

These variations of Tf 1 with 9 and </> help to explain 
the observation that (3 < 1 throughout the entire range of 
T for the polycrystalline sample (Fig. ||). Below ~ 8 K, 
the variation of Tf 1 as a function of (f> causes a very broad 
distribution of Tf 1 in the polycrystalline sample that is 
qualitatively consistent with the established small values 
of p. Similarly, above Tn, the narrower distribution of 
Tf 1 caused by its variation with 9 (Fig. [To]) results in 
a value of (3 that is slightly less than one. In the range 
10-20 K, the moderate range of 9 near 90° selected by the 
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rf pulse at the center of the spectrum and the variation 
of Tn in the sample are probably the major conditions 
responsible for the measured value of (3. 

Our Tf 1 measurements at the highest Bq of 41.5 T for 
a limited number of angles and temperatures are shown 
in Fig. |lj and its inset. No significant deviation from 
the results at lower Bq is evident. This behavior shows 
that although the dynamical properties of the electron 
moments that form the AF state are quite sensitive to 
the alignment of B , they are nearly independent of its 
magnitude up to 41.5 T. 
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FIG. 15: Spin-echo decay as a function of r in a LiVGe20e 
single crystal at 9.0 Tesla. Inset: FFT of the spin-echo decay. 



Figure 15 shows the oscillatory behavior of the ampli- 
tude of the spin-echo decay as a function of the pulse 
spacing r at 9.0 T, 36 K, and 9 = in a single crystal 
sample. We attribute the modulation of the echo height 
to the static quadrupole interaction of the 7 Li nuclei. 
From the period of the modulation, r m , the quadrupole 
frequency vq, is obtained^ vq — l/r m . Since the orien- 
tation of the electric field gradient (EFG) tensor at these 
sites is not known, our limited measurements do not pro- 
vide an exact value of vq. However, from the fourier 
transform of the decay curve for a Ti measurement with 
one of the shortest periods (see inset FigJlJ), an approxi- 
mate value of vq w 15 kHz is obtained for the quadrupole 
frequency. Actually two frequencies vq and 2i/©-are seen, 
as expected for a nuclear spin 1 = 3/2 system.!!!! Also no 
change in the modulation was observed over the entire 
range of T that was covered. This behavior indicates that 
the EFG is constant and that no structural change occurs 
at Tn, in agreement with neutron and x-ray diffraction 
measurements J3 

In Fig. r 2 _1 as function T is plotted at B = 9.0 T 
for both a polycrystalline powder sample and a sin- 
gle crystal sample. Several features are seen in these 
measurements. Except for very close to the transition, 
T 2 _1 ~ 1000 s _1 above Tn- Also, there is a narrow peak 
in T% within 1 K of Tn in which T 2 _1 is enhanced by 
about 25 % for the powder sample and by 70 % for the 
single crystal one. We attribute these increases to slow, 
critical fluctuations, of the local field near Tn- 

As T decreases below 23 K, there is an increase in T^ 1 
by the factor 2.3. From elementary considerations, one 
might expect a reduction in T^ 1 caused by the AF field 
"detuning" the 7 Li spins. Although it should be present, 
this mechanism is clearly not the dominant one for T% ■ 
There are three mechanisms that are often responsible 
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FIG. 16: 7 Li T 2 _1 as function of T at 9.0 T in polycrystalline 
LiVGe2 06 (filled circles) and a LiVGe2 06 single crystal (open 
circles) . 



for such an increase in 2TT 1 : (1) A slow fluctuation of 
the local magnetic field or EFG in the frequency range 
near to T 2 7 , (2) an increase in T^ 1 to values on the or- 
der of or larger than T^ 1 caused by other mechanisms, 
or (3) an enhancement of the spin-spin interactions be- 
tween the nuclei being measured. We expect that the 
first is not determining T^ 1 because it seems unlikely 
that slow fluctuations would be independent of T when 
the fast ones that determine T^ 1 vary so rapidly with T. 
The second obviously does not apply because all of the 
measured values of TT' 1 are too small. It may be that the 
third mechanism does apply through the 7 Li- 7 Li spin- 
spin interaction being enhanced by the Suhl-Nakamura 
interaction mediated by the AF spin wave modeslla that 
are expected to form at low T. Although this approach 
shows some promise, to evaluate it in detail is beyond 
the scope of this paper. 



IV. DISCUSSION 

In this Section, we discuss and interpret features of our 
results that have not been covered in the above presen- 
tation of the data. 



A. Magnetic structure 

In this Section we discuss how the NMR spectra of 
both the powdered and the single crystal samples can be 
used to infer the spatial arrangement of the moments in 
the AF phase of LiVGe20g. First, we consider what can 
be inferred from the data in the randomly oriented poly- 
crystalline sample. Its NMR spectra in the AF phase 
can be analyzed in a similar way as in the paramagnetic 
phase (see Sec. Ill A); i.e. sum the contribution of all 



the V 3+ ions in a sphere of about 5nm around the 7 Li 
nuclei, but with the magnetic moments on the V 3+ ions 
having an AF configuration. We used the AF structure 
corresponding to a simple AF period along the chains, 
a ferromagnetic order between the chains and a value of 
1.14 /xb fof the magnetic moments reported by Lums- 
den et al. ,u but with several different orientations of the 
magnetic moments. For the moments pointing along the 
(100) direction, a total field of only 0.022 T, which is 
much less than the measured value of about 0.106 Tesla 



(see Sec. Ill A), is obtained. However, when the moments 
are parallel to the (001) direction, the result is 0.104 T, 
which is very close to the measured value. Thus, the poly- 
crystal measurements are compatible with the proposed 
AF structure by Lumsden et alia, but with the magnetic 
moments pointing along the crystallographic c-direction 
instead of the a-direction. 
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FIG. 17: Shift of the 7 Li NMR-absorption peaks in 
LiVGe 2 6 as function of 6 at 9.0, 23.0 and 41.5 T. The solid 
lines are fits to the data at 9.0 Tesla (see text) and the dashed 
lines are guides to the eye. Inset: Frequency difference Av 
between the two NMR lines in the AF phase, as function of 
Bo at T ~ 10 K. The solid line is a guide to the eye. 

The measurements on the single crystal samples pro- 
vide even more direct evidence that the hyperfine fields at 
the Li sites, and therefore the moments on the V 3+ ions, 
are aligned parallel (and antiparallel) along the crystal- 
lographic c-direction. This can be seen from Fig. [l7| , 
where Av for the two peaks in the 7 Li spectrum caused 
by the magnetically inequivalent sites in the AF phase 
near 10 K is plotted as function of 9. For the measure- 
ments at 9.0 T, a good fit to Av is given by: 

Av = 149.01 ± 1.86 cos(0) + 0.01 cos(26>) T 0.03 cos(30). 

The dominant term in this fit is proportional to cos(0), 
which is consistent with Bj being parallel to the c-axis 
for all values of 9. Since B; has this orientation for the 
AF ordered moments parallel to the c-axis, we conclude 
that the moments themselves are aligned with the c-axis 
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in the AF phase. 

For the measurements at 9.0 T, the deviations from the 
cosine function are quite small. The cos(26>) dependence 
can ba-interpreted as a slight tipping of the moments 
by Bo-Ej The inset of Fig. [IT] shows the maximum shift 
Az/ max as a function of Bq at T m 10 K. The few high 
field measurements near 10 K indicate a small reduction 
of the order parameter. Because of the small number of 
measurements, this point needs confirmation by a more 
complete set of measurements. 



B. Phase transition 



As shown in Sec. [II A , the NMR spectra of the poly- 
crystalline sample indicate the coexistence of the para- 
magnetic and AF phases over a range of 1.5 K, centered 
around Tn- It has been suggesteda that a first order 
transition is responsible for the coexistence of the two 
phases. On the other hand the T-dependence of the split- 
ting, which is proportional to the order parameter, varies 
continuously and smoothly to zero at the transition tem- 
perature. This behavior indicates that the phase tran- 
sition is of second order, or at most, very weakly first 
order. An alternative explanation for the coexistence of 
the two phases is a distribution of transition tempera- 
tures Tjq. Such a distribution, caused by dislocations, 
stacking faults and V vacancies, could easily be present 
in a polycrystalline sample. This explanation is also con- 
sistent with the small value of the exponent near the 
transition (inset of Fig. ||), which indicates a very broad 
distribution of T^ 1 about Tn, where T^ 1 varies more 
rapidly with T than at temperatures nearby. The vari- 
ation of Tj -1 with T very close to the transition (inset, 
Fig. ||) is slower than expected given the divergent criti- 
cal behavior near a second order transition. As indicated 
earlier, this suppression of the critical divergence could 
be caused by a distribution of Tn. 

This interpretation in terms of a distribution of Tn 
is supported by our measurements on single crystals. In 
comparison with the polycrystalline sample, the data im- 
ply a coexistence of both phases over the substantially 
narrower range of about 0.3 K and a peak in T^ 1 that 
is narrower and higher near Tn (inset, Fig. [To| ). This 
behavior is qualitatively consistent with a distribution of 
Tn in which the single crystal has fewer imperfections 
and therefore a more narrow distribution of Tn than the 
polycrystal line s ample. 
Sec. " 



In 



IV A 



it was shown that the AF state of 
LiVGe206 has a rather simple magnetic structure. It 
has, however, some unusual features which we turn to 
now. First, consider the magnitude of Tn. We start with 
a simple spin Hamiltonian which, as will be discussed 
later might not be sufficient to describe the system: 

h = J,, s i ■ s i + D E( 5 D 2 + j ±J2 Si - s i > ( 8 ) 



and denotes an intrachain one. Jn and J± are re- 
spectively the the intrachain and interchain coupling con- 
stants, and D is the single-ion anisotropw. The value 
about 45 K has been estimated for J||.Ell3 Because of 
the quasi 1-D character of the crystal structure, it is 
surprising that the measured value Tn ~ 25 K could 
be more than half J||. On the basis of a mean-field 
calculations and the assumptions Jii = 45 K, four neigh- 
boring chains, and Tn = 25 K, J± has been estimated! 
to be aboutl-4 K, or — 0.03. More advanced 

calculationsEJiijO show, however, that such a ratio is 
barely enough to induce an AF transition, and very un- 
likely to have a value of Tn as large as Jii/2, In fact 
based on the calculation of Senechal et al.t£l we esti- 
mate the ratio of the interchain to the intrachain cou- 
pling constant which would be needed to explain the rela- 
tively high transition temperature to be larger than 0.06. 
From crystal considerations, however, this value seems to 
be very large. For example, AgVP2Sg, which is a com- 
pound with a steucture similar to LiVGe2 06, has a ratio 
Jj-/J\\ < 10~ 5 Jla Therefore, it seems that the relatively 
high transition temperature cannot be explained by the 
interchain coupling alone. If, however, the Haldane gap 
is not present, even a very small coupling between the 
chains might be enough to induce an AF phase transi- 
tion. Later, we will consider this possibility in terms of 
a large single-ion anisotropy. 

Another remarkable feature of LiVGe206 is the weak 
influence of Bq on the properties of the AF state. This 
applies to Tn, is seen in the properties of an apparent 
gap in the fluctuations responsible for (to be dis- 
cussed later), and is evident in the absence of a spin- flop 
transition for Bo as high as 44.5 T. The last point is 
noteworthy because on the basis of a simple mean-field 
apjproximation, the spin-flop field (Hsf) is expected to 



Hsf — (2#e-Ha — H 



2x1/2 



(9) 



where denotes an intrachain nearest-neighbor pair 



where Ha is the anisotropy field and Hf, is the exchange 
field. The field Ha is related to the single-ion anisotropy 
according to D — Ha9^bI 'S. From these considerations, 
D > 50 K, aihich might be responsible for closing the 
Haldane gapEJ and may explain the relatively high AF- 
phase transition temperature Tn- Although the mean- 
field approximation used here might not be fully appro- 
priate because Ha is comparable to Hf, it does allow 
us to obtain at least a rough estimate of the single-ion 
anisotropy using the Hamiltonian of Eq. ^. 

It should n be pointed out that according to a recent 
publication,!] no phase transition is expected to occur in 
LiVGe206- In this work, mid-gap states are assumed 
to be responsible for the susceptibility anomaly in the 
experimental data of LiVGe206 and it is predicted that 
this anomaly will be weaker if their are fewer crystal de- 
fects and non-magnetic impurities. This interpretation is 
in clear contradiction with the prior NMR0 and neutron 
diffractioncl experiments and our NMR measurements re- 
ported here. Our experiments on both polycrystalline 
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and high quality single crystal samples show clearly that 
a magnetic phase transition occurs and permit us to re- 
fine the magnetic structure proposed earlier 3 Although 
the phase transition appears to be more complicated than 
a simple ordering driven by inter- and intrachain cou- 
pling, there is no doubt that the V moments order anti- 
ferromagnetically below approximately 25 K. 



C. Tf 1 as function of T,B, 9 and <f> 

Although in principle part of the coupling responsible 
for Tf 1 of the 7 Li spins could be quadrupolar, the dis- 
cussion in the next paragraph argues that there is no sig- 
nificant quadrupolar contribution to it. As a result, our 
interpretation of Tj -1 involves only magnetic coupling to 
the Li spins. 

There are basically two mechanisms for electric field 
gradient fluctuations that might contribute to T^ 1 of 
the 7 Li spins: Charge fluctuations associated with some 
instability of the lattice and the Raman-type phonon 
process first described by van KranendonkBj Since the 
quadrupolar splitting found in Sec. IIIE| is very small 
and does not change over the temperature range stud- 
ied, there is no reason to expect a significant contribu- 
tion to T^ 1 from charge fluctuations. Also, the small 7 Li 
quadrupolar interaction and observed T-dependej*6e of 
T^ 1 rule out quadrupolar relaxation by phononsEil We 
therefore conclude that there is no significant quadrupo- 
lar contribution to T^ 1 and that it is caused by magnetic 
fluctuations. 
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FIG. 18: 7 Li Tf 1 as function 1/T in polycrystalline 
LiVGe206 at 9.0 and 23 T. The solid dashed lines are fits 
to the low T data using Eq. [Io| The values shown for A are 
14 K at 9 T and and 11 K at 23 T. 

In Fig. [l^thc spin-lattice relaxation rate, T^ 1 , of poly- 
crystalline LiVGe2 0@ is plotted as function of T _1 for 
B Q = 9.0 T and 23 T. Well into the AF phase below 10 K, 
the behavior of T^ 1 has the T-dependence expected of 



electron spin excitations across an energy gap (A); i.e. 
the slope of the curve is constant and negative at low T. 
If one simply fits the data to 



T^\T) cx exp(-A/T) 



(10) 



the values obtained for A are 14 K at 9.0 T and 11 K at 
23 T. 

Such values must, however, be interpreted as an aver- 
age over a distribution of A that is quite broad. This is 
seen in Fig [l^, where there is a large variation of the slope 
for two different values of <j> at 28.5 T (solid lines), and 
a large slope at 9 T and ~ (dashed). The two values 
A =2.6 K and A =14.7 K indicate a very large effect 
of the alignment of Bq in the azimuthal plane on A for 
magnetic excitations at low T. Therefore, the following 
discussion of the polycrystalline sample results applies to 
an average behavior and is approximate and phcnomeno- 
logical. 

First, consider the effect of Bq on A. In a 3-d antiferro- 
magnet, when there is a gap in the excitation spectrum, 
it usually depends strongly on B because of the Zeeman 
interaction, which is oc q^lbBq. Instead, our measure- 
ments show the very weak field dependence of less than 
3 K for a difference in applied field of 14 Tesla. Since 
g(i B B/k B , with g = 1.79, is about 28 K for B = 23 T, 
the usual expectations for the dependence of A on Bq do 
not apply. 

A ap- 
From 



Also, in comparison to our results in Sec. IVB. 
pears too small to be attributed to a magnon gap. 
A = 14 K and 11 K at 9.0 T and 23 T respectively, one 
would expect a A ~ 16 K. From this value, a spin-flop 
field smaller than 20 T is expected, which is contradicted 
by the absence of a spin-flop transition up to 44.5 T 
shown by our NMR-spectrum measurements. On the ba- 
sis of this evidence, it appears that at low T the magnetic 
fluctuations responsible for T± in LiVGe2C>6 are not 
simple magnon excitations. This is particularly evident 
from the Tj -1 measurements shown in Figs. [H] and 
where very large anisotropics in T^ 1 are seen. 

It is rather difficult to identify the microscopic mech- 
anism responsible for the very large anisotropy of T^ 1 
shown in Fig. [l^, where, for comparison, the variation 
cos 4 4> and cos 8 </> are shown by the dashed and solid lines, 
respectively. For example, if it were caused by a fluctu- 
ating magnetic field aligned with the value of <j> corre- 
sponding to the maximum in , for other values of 
one would expect the much weaker angular variation 
cx cos 2 4>. Similar arguments for quadrupolar relax- 
ation (excluded above from other arguments) by a fluc- 
tuating EFG could give a variation up to cos 4 (f>. 

Although we do not have a microscopic model for 
it, the temperature and angular variations seen in 
Figs. |l2| and |l3] suggest a gap- type behavior with a split- 
ting which itself is very anisotropic. More work, both 
experimental and theoretical, is needed to identify the 
mechanisms for this behavior. 

In summary, the angular dependence of the spin-lattice 
relaxation rate measurements for the single crystal sam- 



Tf 1 
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pies is rather complicated and it is difficult to construct 
a detailed interpretation. For both the paramagnetic 
regime and the AF regime down to about 8 K, T^ 1 
Fig. [ll] and Eq. | indicate that the largest contribution to 
T^ 1 has the angular dependence oc sin 2 9. This behavior 
is consistent with magnetic fluctuations that are predom- 
inantly along the c-direction that have dipolar coupling 
and only a small isotropic contribution. It may reflect 
primarily amplitude fluctuations of the AF order param- 
eter. 

Below about 8 K, the disappearance of the 9- 
dependence of T-j -1 and the emergence of its in- 
dependence indicates that at low T, the origin of the 
fluctuations responsible for T^ 1 is very different from 
what it is at higher T. 



D. Orbital degrees of freedom and magnetic 
anisotropy 

In this section we mention some points about the mag- 
netic anisotropy and the orbital degrees of freedom of 
the V electrons in LiVGe20e suggested by our measure- 
ments. They suggest that the second order splitting 
Acf of the t2 g orbitals is rather small, with the con- 
sequence that (i) there is a large uniaxial anisotropy as 
<7_i_ = 2(1 — A/Acf) (A = spin-orbit coupling) is strongly 
reduced for small values of Acf, and (ii) orbital fluc- 
tuations may play a significant role in the properties of 

IT 1 . 

A large anisotropy energy has not yet been observed di- 
rectly from 51 V NMR measurements because we were not 
able to observe the signal. It can, however, be inferred 
from several aspects of our 7 Li NMR measurements. In 
particular, the absence of a significant £?o-dependence of 
Tn (Fig. [?]), the absence of a spin- flop transition for Bo 
up to 44.5 T, and the 9 dependence of Au (Fig. |l7|) all 
indicate a large uniaxial anisotropy for the static magne- 
tization. A similar picture emerges from the mea- 
surements above about 8 K, where the ^-dependence of 
T^ 1 (Fig. [ll]) indicates that the fluctuations of the V 
moments are constrained mainly to the c-direction. 

The situation of the V ions in LiVGe2Q« is similar to 
those in V2O3, where it has been reportedE] that, at least 
for the metallic phase, T^ 1 at the V sites is dominated by 
orbital fluctuations. In LiVGe206 the presence of orbital 
fluctuations of the V 3+ t 2g orbitals is suggested by the az- 
imuthal 4> dependence of the spin-lattice relaxation rate 
at the 7 Li site (Figs. [l2| and |l3|) which is very anisotropic 
and independent of the magnitude of Bq. Furthermore, 
orbital fluctuations might be very effective at the 51 V site 
itself, and may be responsible for the absence of the 51 V 
NMR signal. 

According— to a quantum chemistry analysis of 
LiVGe20 6 ,Era if A CF between the V 3+ orbitals is simi- 
lar to the hopping integrals between neighboring V sites, 



it is not possible to describe the system with a pure spin 
Hamiltonian and the orbital degrees of freedom have to 
be included explicitly. Our NMR data indicate that this 
is the case in LiVGe2 0g. Therefore, the Hamiltonian 
given in Eq. ^ should be modified to better describe the 
physics of LiVGe2 06, and include orbital degrees of free- 
dom which might or might not be strongly coupled to the 
spin degrees of freedom. 



V. CONCLUSIONS 

We have presented pulsed 7 Li NMR measurements in 
polycrystalline and single crystal samples of the quasi 
one-dimensional antiferromagnet LiVGe2 06 over the Bq 
and T ranges 9-44.5 T and 1.5-300 K respectively. They 
cover both the paramagnetic and the AF phases, for 
which the transition is at Tn — 24.5 K. The measure- 
ments include NMR spectra and the relaxation times Tj -1 
and T 2 -1 , often as a function of alignment of Bq. From 
the spectrum measurements, we find that in the AF phase 
the magnetic structure is consistent with that reported 
on the basis of neutron diffraction measurements, B but 
with the moments aligned parallel to the c-axis. Mea- 
surements of T 2 _1 show oscillations caused by the static 
electric field gradient. The corresponding interaction is 
quite small and independent of T, which indicates that 
over the range of T that was covered, no lattice struc- 
tural transition is observed in LiVGe2 0g. The spectrum 
measurements also provide the T-dependence of the or- 
der parameter and show that the transition is either sec- 
ond order or weakly first order. The coexistence of the 
two phases over a narrow range around Tn and the be- 
havior of the NMR linewidth below it is attributed to a 
distribution of Tn in the samples. Both the spectra and 
the angular dependence of indicate that the external 
field has at most a small effect on the alignment of the AF 
moment. There is no spin-flop transition up to 44.5 T. 
These features show that there is a very large anisotropy 
energy in this material and that the Hamiltonian should 
include orbital degrees of freedom to adequately describe 
it. Below 8 K, a rapid dependence of T-f 1 on the az- 
imuthal angle suggests the presence of a low energy gap 
for magnetic fluctuations that is highly anisotropic. 
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